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ABSTRACT 

An atomization model prescribing jet breakup when jet distortion caused by aero- 
dynamic forces attains a critical value is used to analytically and numerically evaluate 
the atomization rate perturbations caused by acoustic oscillations. The varied dynamic 
response within traveling, standing, and radial transverse acoustic modes is presented. 
An amplifying effect of a steady tangential velocity on dynamic response is also examined 
and compared with experimental evidence. The results generally imply that jet atomiz- 
ation dynamics can significantly affect the stability of a rocket engine combustor. 
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ANALYSIS OF THE DYNAMIC RESPONSE OF LIQUID JET 
ATOMIZATION TO ACOUSTIC OSCILLATIONS 
by Marcus F. Heidmann and John F. Groeneweg 
Lewis Research Center 

SUMMARY 

The dynamic behavior of an axial liquid jet within traveling, standing, and radial 
transverse acoustic modes in cylindrical chambers is examined analytically and numeri- 
cally for dynamic properties useful in characterizing a combustion process for stability 
analysis. An atomization process is used which prescribes jet breakup when jet distor- 
tion caused by aerodynamic forces attains a critical value. Two response factors are 
evaluated: the in-phase Nj^ and out-of -phase Nj components of the atomization rate 
perturbations with respect to pressure perturbations. Effects of acoustic oscillations on 
jet breakup time are also determined. In general, the results imply that jet atomization 
dynamics can significantly affect the stability of a rocket-engine combustor. 

The response factors Np and Nj are converging periodic functions of a charac- 
teristic breakup time wr with peak value dependent on acoustic conditions. Consid- 
ering pressure oscillations alone, peak values approach unity. Including a velocity per- 
turbation in-phase with the pressure elevates the peak values by a factor of 3/2; an out- 
of -phase velocity perturbation reduces the peak by 1/2. Velocity oscillations also re- 
duce the mean jet breakup time and, thereby, affect the characteristic time off of the 
process. 

Spatial variations of acoustic properties within a given mode can cause large vari- 
ations in local dynamic properties of the atomization process. The effect on the re- 
sponse factors of the varied velocity and pressure patterns within a specific mode pre- 
cludes any simple characterization of dynamic properties. However, the response fac- 
tors for the traveling mode can be larger than those for either the standing or radial 
modes. 

The amplifying effect of a steady vortex velocity on the response in a traveling mode 
is substantial. Peak values of the response factor can increase by an order of magni- 
tude. The results appear to be consistent with experimental studies of combustion in- 
stability with such vortex flow. 



INTRODUCTION 


In the theoretical analysis of acoustic mode instability in rocket-engine combustors, 
the dynamic properties of the combustion process represent important initial assump- 
tions which substantially affect the solutions for system stability. Also, in the phenom- 
enological approach to instability, the combustion process represents the energy source 
which sustains the acoustic oscillations and becomes a prime target in any attempt to 
stabilize the system. Characterization of the dynamic behavior of the combustion proc- 
ess, therefore, is not only useful but at times must precede an attack on the problems of 
instability. 

Complete characterization of the many facets of the actual combustion process in a 
rocket engine and the integration of them into a comprehensive model of the overall sys- 
tem appears improbable. Instead, dynamic behavior of one or a combination of several 
processes comprising the overall combustion process is frequently employed to simulate 
combustion dynamics. Atomization, vaporation, mixing, chemical kinetics, drop shat- 
tering, and recirculation are among the individual processes which require attention in 
liquid-propellant engines. Theoretical analyses of the dynamic behavior of individual 
processes are particularly helpful in deciding which process or combination of processes 
can best characterize combustion. The analysis of droplet vaporization (refs. 1 and 2), 
for example, has shown a significant response of the vaporization rate to acoustic oscil- 
lations. This response, coupled with the concept of a vaporization limited combustion 
process, can be used to characterize the combustion process in certain engines. 

Similarly, the analysis of the response of a gaseous hydrogen injection system to 
acoustic oscillations combined with the oxygen vaporization response (refs. 5 and 6) also 
provide physical insight into the dynamic properties of the combustion process . In this 
report an attempt is made to expand this insight of combustion dynamics by analyzing 
the dynamic behavior of the breakup or atomization of liquid jets. 

Perturbations in atomization rate caused by acoustic oscillations are of particular 
importance when, as observed in some experimental tests (refs. 7 and 8), extremely 
small drops are formed which burn rapidly. The atomization process may then be the 
rate controlling process in the dynamic system. The low surface tension of such liquids 
as oxygen or severe pressure oscillations (produced naturally or by explosives used in 
dynamic testing) are two conditions conducive to the formation of small drops and com- 
bustion instability dependent on atomization dynamics. 

The liquid jet breakup or atomization model developed by Clark (ref. 9) is used for 
the analysis of jet dynamics. Although elementary in concept, the model expresses the 
instantaneous jet properties during breakup - a necessary prerequisite for an unsteady 
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analysis of any process. The analysis emphasizes transverse acoustic modes in cylin- 
drical chambers because of the prevalence of such modes in unstable rocket combustors. 
Approximate analytical solutions of the response of atomization rate to pressure oscilla- 
tions are derived. These results are supported by numerical solutions which eliminate 
the simplifying assumptions used in the analytical procedures. The effects of a spatially 
distributed breakup process and vortex flow on dynamic behavior are also examined. 



Hydrogen 


Camera 


Jet orifices 


Injector ring 


Condensing lens 


Light source 


Figure 1. - Two-dimensional circular combustor and photographic arrangement. 
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OBSERVATIONS OF JET DYNAMICS 


A review of some experimental observations of a liquid jet in a two-dimensional cir- 
cular combustor will illustrate the type of dynamic behavior of interest in this analysis. 
Figure 1 typifies the experimental arrangement (ref. 8). Injected propellants and com- 
bustion gases flow radially inward from the circumference, and photographic observa- 
tions are centered on a single injector element. Vortex flow within the cavity, induced 
by the tangential injection of a secondary gas, caused unstable combustion coupled to the 
first traveling transverse acoustic mode. 

Typical examples of the periodic behavior of a liquid oxygen jet burning with gaseous 
hydrogen during unstable combustion (refs. 7 and 8) are shown in figure 2. Sequences 
of high-speed motion picture frames are shown where successive frames are time sep- 
arated by a period differing somewhat from the oscillating period of the instability 
(a stroboscopic effect). Figure 3 illustrates the corresponding acoustic conditions for 
successive frames in a typical sequence. This stroboscopic technique resolves gross, 
rather than detailed, changes in the jet shape with time. Examination of such pictures, 
however, showed a time varying jet length which is coupled to the acoustic oscillations. 
Such jet-length variations imply variations in mass release rate. These observations 
led to the suspicion that unstable combustion was directly related to the dynamic behavior 
of the liquid jet. 

Experiments using mixed oxides of nitrogen and hydrazine as propellants in the two- 
dimensional combustor also implicated liquid jet behavior as vitally important during un- 
stable combustion. In these tests (ref. 10), the impingement of unlike propellant jets 
was used for injection. Although combustor stability characteristics were generally 
similar to the hydrogen-oxygen test, the occurrence of a nonlinear instability is of par- 
ticular interest. Although the general behavior was a linear instability with pressure 
amplitude related to the degree of vortex flow, a nonlinear instability independent of the 
vortex flow was triggered above a critical pressure amplitude. The liquid jets were 
photographically observed in the neighborhood of this threshold region separating linear 
and nonlinear instability. The observations showed that nonlinear instability occurred 
when the acoustic oscillations were of sufficient magnitude to cause jet breakup prior to 
the impingement point. The stroboscopic sequence of photographs in figure 4 show a 
condition of nonlinear instability. In this case, breakup of the outer (mixed oxides of 
nitrogen) jets prior to the point of impingement appears to be a critical condition. The 
severity of this nonlinear instability was reduced by protecting the jets with thin walled 
tubing and thus maintaining impingement. These tests imply that dynamic behavior of 
the jet breakup process determines not only the threshold condition for nonlinear insta- 
bility but that the acoustic gains associated with simple jet atomization are larger than 
those associated with impinging jet atomization. 
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(a) Pressure amplitude, 20 percent; low oxygen jet velocity. 



(b) Pressure amplitude, 10 to 20 percent. 



(c) Pressure amplitude, 5 to 10 percent. 



(d) Pressure amplitude, 5 percent; high oxygen jet velocity. 

Figure 2. - Liquid oxygen jet behavior with first transverse mode of instability (ref. 7). 
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Figure 3. - Velocity vector and phase position for sequence of jet photographs (ref. 8). 
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(a) Injection, fuel-oxidant-fuel; stable combustion. -H N- 0.1 in. (2.54 mm) 



(b) Injection, fuel-oxidant-fuel; self-coupled oscillation 



(c) Injection, oxidant-fuel-oxidant; stable combustion. 



(d) Injection, oxidant-fuel-oxidant; self-coupled oscillation. 

Figure 4. - Impinging jet behavior with first transverse mode instability. Mixed oxides of nitrogen and hydrazine propellant combination. 
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The tests with the two-dimensional combustor emphasize the need for an analytical 
description of liquid jet dynamics. The experimental tests provide only qualitative de- 
scriptions with perturbations in mass release rate, the crucial quantity for instability, 
being most elusive. The physical insight gained from such tests, however, is valuable 
in formulating an analytical model of jet dynamics. 


LIQUID JET ATOMIZATION MODEL 

In the tests with the two-dimensional circular combustor, the transverse acoustic 
particle velocity exhibited a dominant effect on the breakup of the liquid jet. The effect 
of a transverse flow of gas on jet breakup, therefore, is of particular concern in any 
dynamic analysis. 

The effect of a steady transverse flow of gas on the breakup of water jets was studied 
by Clark (ref. 9) . In that study, an analytical model was formulated to support the ex- 
perimental effect of steady flow on jet breakup. The model is particularly useful in ex- 
amining the dynamic behavior of a jet. It conveniently describes the instantaneous jet 
properties during breakup and is not restricted to integral properties such as jet length 
or breakup time as are some other models (which restrict their dynamic implications). 

Figure 5 illustrates the model formulated in reference 9 for the basic mechanism of 
breakup of a liquid jet element in a cross flow of gas. The mechanism is similar to that 
used and observed for drop breakup in reference 11. 

The external pressure distribution creates a pressure gradient within a element of 
liquid jet which flattens the liquid element. At the edges of this flattened element drops 
and ligaments are torn off by the combined action of the tangential stress and surface 
tension. Breakup of this type is restricted to conditions where external forces are large 
compared with surface forces. 





Undistorted 
cross section 


Distorted cross Breakup 

section 


Figure 5. - Basic mechanism of liquid jet element breakup (ref. 9). 
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The distortion 6 of an element of jet caused by an acceleration a persisting for a 
period t from the time of injection t - r is given by 


6 = 



( 1 ) 


The acceleration is caused by the external forces acting on the element of liquid 
mass. From a consideration of the pressure distribution around the cylindrical element, 
the acceleration becomes 


Pi D 


( 2 ) 


Breakup occurs when the distortion relative to the element diameter e = 6/D attains 
some critical value. The basic mechanism for breakup, therefore, is defined by 


=-— r f Pgv 2 dt ; 

-.2 •'t-T Jt-T £ 

P;D 


** dt x 


(3) 


For a steady flow of gas the breakup time r__ given by the integration of equa- 

ss 


tion (3) is 


ss V W 


,V2 


(4) 


Steady flow experiments with water jets showed that an e of 3 to 4. 5 characterized a 
mean breakup time r gg where only 50 percent of the jet element mass remains intact. 
Complete breakup was characterized by an e of 10 to 15. 

The study of Clark revealed other properties that are useful in the dynamic analysis 
of jet breakup. Although equation (3) describes the basic mechanism, the actual breakup 
process exhibited distributed properties that are not implied by the equation. In particu- 
lar, the breakup of a jet did not occur after a discrete time t from injection as implied 
by equation (3). Rather, breakup or mass removal occurred over some finite time or 
along some finite length of jet. On the average, the measured distribution of the mass 
release rate as determined from figure 17 of reference 9 was distributed with time or 
length as shown in figure 6. These distributed properties will be considered in the dy- 
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namic analysis. Initially, however, the discrete breakup mechanism given by equa- 
tion (3) will be examined for dominant dynamic properties. The effect of the distributed 
nature of breakup on dynamics will be examined separately. 

The study by Clark also showed that the velocity V should be the gas velocity rela- 
tive to the liquid element. Expanding this concept to an axial jet in a cylinder, the mag- 
nitude of the relative velocity vector acting on an element of jet is given by 

V= (l u i" u Z | 2+u 0 + u r) / 

This definition of the effective velocity changes the concept of simple element flattening 
previously conveyed. Rather, jet distortion involves a more complex motion of liquid 
mass caused by cummulative momentum transfer from the gases independent of the di- 
rection of the applied force. Clark’s experiments, although not complete, verified a 
breakup time independent of force direction. Proceeding with this concept, the basic 
mechanism for jet distortion caused by transverse acoustic oscillations, where u^, u r , 
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and p 

S 


are time varying quantities, becomes 


£ 



rt rt x 

f / f(t xx )dt xx dt x 
A-t A-r 


( 6 ) 


where 


f(t) = p g (t) []u z - u z | 2 + u 2 (t) + u 2 (t)] (7) 

Equation (6) is the analytical expression of the breakup mechanism used herein to 
examine the dynamic properties of liquid jet atomization. Dynamic properties are pri- 
marily established by the time dependent aspects of equation (6) . For comparison, the 
approach of Buffum and Williams (ref. 12) for the analysis of jet deflection would give an 
expression with similar time dependent properties when breakup is characterized by a 
critical value of jet deflection. 


THEORETICAL ANALYSIS 

Analytical expressions based on certain simplifying assumptions required for closed- 
form solutions will be derived for the dynamic response of jet atomization in a trans- 
verse acoustic mode. 

The atomization rate w’ for a jet with a discrete breakup time t and a constant 
velocity is given by 


w*(t) = - HI (8) 

dt 

An expression for dr/dt for jet breakup in an acoustic field is obtained from a con- 
sideration of equation (6). In a time varying environment, each element of jet experi- 
ences a different time history of applied forces . If the criterion of a nonvarying critical 
distortion e for all elements is applied, the following must be satisfied 


d£=0 = -2- /■* r’Vldt-dt* (9) 

dt Pj D 2 Jl - T Jt ~ T 
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Applying the Leibnitz rule for differentiation of an integral results in this expression 
of equation (9): 

J rt 

' f(t)dt =0 (10) 

t-T 

Reapplying the Leibnitz rule to the first term of equation (10) gives the expression 

-(l - — V(t - t) dt + f * f(t)dt = 0 (11) 

\ dt/ Jt-T •'t-T 

or 

-r|l ” d^ f ^ ~ ^ + Jt = ° U 2 ) 

Solving for -dr/dt gives the atomization rate w’ as 



w 

The function f(t) depends on specific properties of an acoustic mode at the location 
of the jet. Equation (7) uses a conventional expression for the velocity components. For 
convenience in analyzing the general case of transverse acoustic modes, however, a dif- 
ferent notation will be used. The transverse velocity will be resolved into two phase- 
related components. These are the in-phase u^ and the 90° out-of-phase Uj compo- 
nents with respect to the time variations in pressure. These components can be readily 
identified from the usual first-order solutions for acoustic modes. In this notation the 
function f(t) is given by 

f(t) =p g (|u z - u z | 2 + U R+ u i) (14) 

For any of the transverse modes, the pressure, density, and velocities are of the 
form 


i _ _ dT _ 

dt 


r 

Jt-T 


f(t)dt 


rf(t - r) 


- 1 


(13) 
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(15) 


P(t) = P(1 + P r cos cut) 


p g (t) = 


= Pg — P' cos cut) 


u R (t) = u R cos cut 


Uj(t) = Uj sin .cut 


(16) 

(17) 

(18) 


Using these properties, the function f(t) expressed as a harmonic series becomes 
i(t) = Pg [(l V 1 2 + \ + \ S I 2 ) + 1 S' (l V \ 1 2 + | Sr + ~ 50COS «t 

+ 2 (^H ” ^l) cos 2wt + ~ P’(u R - Uj 2 jcos 3cutj (19) 


Neglecting the higher harmonic terms, the function takes the form 


*<*> = Vrmsl 


1 + V’ — cos cut] 

v y ) 


where 


V’ = 


I 1 2 3-2 1-2 

h- u J + j u r + j u i 

|u, -u | 2 +ic| + IS? 
1 1 z 1 2 R 2 1 


( 20 ) 


( 21 ) 


and 


, r 2 | 1 2 1-2 1-2 
V rms= K’V + 2 U R + i U I 


( 22 ) 


Performing the integration indicated in equation (13) with equation (20) as the forcing 
function results in the solution obtained by assuming that the difference between r and 
r is small: 
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(23) 


** T>T 

w'(t) = V — 


sin cot — 4; sin w (t - t) - cos co(t - f) 


(jOT 


(jOT 


1 + V* — cos co(t - t) 


y 


Figure 7 shows several examples of the atomization rate perturbations given by 
equation (23) for several values of wr. Although the atomization rate oscillations are 
not pure harmonic functions, they show a change in phase angle and amplitude relation 
between pressure and atomization rate with a change in breakup time. 



Figure 7. - Effect of mean breakup time on atomization 
rate perturbations for amplitude parameter, 

(V 1 / y)P\ of 0. 1. 


One method of examining such dynamic behavior is to extract two response factors: 
the in-phase N R and out-of-phase Nj components of the atomization rate with respect 
to pressure. Since the pressure oscillation was assumed a cosine function, the first- 
order terms for the atomization rate are given by 


w’(t) = I”(N r cos cot + Nj sin cot) 


(24) 



where the response factors are given by 


and 



1 

fr(P’) 2 



w’(t)P' cos cot dcot 




w'(t)P' sin wt dwt 


An alternate expression for the first-order expression is 

w'(t) = w' cos cu(t + 9) 
where the gain and angle are defined by 


(25) 


(26) 


(27) 


w 


P’ 



> 


9 = tan 


-1 



(28) 


In either of these notations, and Nj represent the real and imaginary part of a 
type of admittance relating atomization rate and pressure. 

Using the atomization rate given by equation (23), the integral equation (25) and (26) 
define and Nj for jet atomization in a transverse acoustic mode. Although no exact 
analytical solution exists for these integrals, an expansion to first order in P’ gives the 
result, 


and 



(29) 


(30) 
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(b) Gain and phase angle. 

Figure 8. - Dynamic properties of jet atomization pro- 
cess from theoretical analysis. 


Figure 8 shows the dynamic properties of the atomization process given by these 
solutions. Both and Nj are periodic with off and exhibit extreme values of about 
±( V' / y) (fig. 8(a)). In general when either Npj or Nj is a maximum or minimum, the 
other is approximately equal to zero. Conditions at c or near zero are an exception. In 
this region, the gain, w'/l?', decreases with c ot (fig. 8(b)) and causes both and Nj 
to approach zero. When off is greater than tt, the gain is relatively constant but the 
phase angle 9 increases nearly linearly with off. 

The response factors and Nj are useful boundary conditions representing the 
combustion process in a stability analysis of a combustor (ref. 13). The quantitative 
significance of these properties on stability can only be established from a system anal- 
ysis that mates combustion properties with system acoustics. Qualitatively, the compo- 
nent of the atomization rate in phase with the pressure oscillation is of deciding im- 
portance in many instances. It is of primary importance in the phenomenological ap- 
proach to combustion instability. If the magnitude of N-^ of a process approaches unity, 
the contribution of the process to unstable operation is substantial. 

The magnitude of for the atomization process can fall in the critical region 
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near unity. The maximum value of N R is about V’/y, where V f is given by equa- 
tion (21). The maximum value of N 0 for the condition of an axial velocity difference 

alone (u R and Uj = 0) is 1/y. For u R alone, (N R j is g for u i alone, 

( f \ j max 

Nr] is ^ y. These values are near unity and indicate that the jet atomization 

process can be an important dynamic component of a combustion system. 

In addition to these effects of acoustic velocities on the magnitude of the response 
factors, the mean breakup time of the jet is also velocity sensitive. This can be shown 
by considering only the first harmonic content of f(t) given by equation (20) and perform- 
ing the integration indicated by equation (6). The limiting solution as c of becomes 
large is 



This form of the solution is identical to that used by Clark for a steady flow (eq. (5)). 

Acoustic affects enter into the definition of the velocity term V„ m _ given by equa- 

i m s 

tion (22). The effect of velocity perturbations on mean breakup time is shown in figure 9. 



0 100 200 300 

Acoustic particle velocity, (Gr + Qj) 1/2 , m/sec 


Figure 9. - Effect of acoustic velocity oscillations on 
mean breakup time. 
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An increase in velocity perturbations decreases the breakup time below its steady-state 
value. The effect decreases with an increase in |u^ - u | . 


NUMERICAL ANALYSIS 


The theoretically derived dynamic properties of jet atomization involved a number of 
simplifying assumptions of unknown quantitative significance. In general, the procedures 
resulted in analytical limitations similar to those restricting linear analysis to small 
perturbations. An important application of jet dynamics, however, is in the region of 
finite and large acoustic oscillations. These uncertainties of the analytical procedures 
prompted a numerical analysis which preserved all nonlinear behavior of the assumed 
mechanism of atomization. An analysis based on the following procedure was pro- 
grammed for a digital computer to quantitatively establish nonlinear dynamic properties 
for some representative acoustic conditions. 

The basic breakup mechanism described by equation (6) and the general description 
of f(t) for a transverse acoustic mode given by equation (19) were assumed. An exact 
integration of equation (6) between the time of injection t^ = t - r and the time t was 
performed for these acoustic properties. The solution which defines the distortion of any 
element of jet is given by 


\ CJ ^ A 

— — u> j -i e = u> 2 (t - t Q ) 2 + 2 — V' [cos wt 0 - cos wt - w(t - t^) sin wt 0 "| 

In v L -* 


rms 


j (u| - uf\ r 

+ — i cos 2wt n - cos 2ont - 2co(t - t n ) sin 2wt n 

4 / \2 L 0 0/ °J 

Y rmsj 

. A 2 _ -2\ 

+ — — L l [cos 3o>t n - cos 3wt - 3o>(t - t n )sin 3a>t n "| 

18 y / y \2 L 0 0 0J 

y v rms ) 


(32) 


where V and V’ are defined as previously in equations (21) and (24). 
rms 

On specifying acoustic properties for a particular mode, tables of values relating 

f(D/V H 2 (p 7 /p )e to wt can be generated for progressive incremental values of 
L rms J t ^ 

cotp encompassing one cycle of oscillation. The breakup criterion specifies that any 
element will breakup when e and, therefore, [(D/V rmg ) cJ] 2 (p^/Pg ) e reaches some con- 
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stant value. On specifying a specific value of this modified distortion parameter, a 
search of the tables will give interpolated values of cut which satisfy the breakup cri- 
terion for progressive values of cutp. The difference cut - cutp is the dimensionless 
breakup time cot for successive jet elements. From such results, the atomization rate 
w'(t) is given by the incremental change -(Acur/Acut) as indicated by equation (8). These 
values of w’(t) used in a numerical integration of equations (25) and (26) give the in- 
phase and out-of-phase Nj response factors. A value of c of characterizing a 
selected value of the modified distortion parameter is established by a numerical average 
of cur for the elements contained within one cycle of oscillation. 

In the actual procedure, tables of modified distortion parameters for time incre- 
ments of 1000 per cycle were computed for 100 values of initial time. 


APPLICATION TO SPECIFIC MODES 

A general description for acoustic properties was used in the development of the 
theoretical and numerical analysis of jet atomization. Application of these analysis to 
several specific conditions of acoustic resonance will illustrate some of the important 
effects of variable acoustic properties on jet dynamics. Three transverse acoustic 
modes will be considered. These are the first or lowest order traveling, standing, and 
radial modes. 


Traveling Transverse Acoustic Modes 

The first-order terms for the properties of a traveling mode are given by 

P = pjl + AyJj(o)cos(cut + 0)j 

P g = P g [l + AJ 1 (o)cos(cut + 0)] 


J i (°0 

u = Ac cos(cut + 9) 

a a 


u r = Ac 


Jn(«) 


J,(<*) 


a 


sin(cut + 9) 


where a = 1. 841 r. 


(33) 

(34) 

(35) 


(36) 
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These properties establish the following identities for the notation used in the gen- 
eral description of a transverse mode in equations (15) to (18). 


u„ = u. 


R 




u = u T 
r I 


P’ = AyJ 1 (Q!) 
CP’ 


U R * 


Oly 


U I = 


CP' 

OLy 


a 


Jq( 

J,(«) 


y 




(37) 


The principal features of this traveling mode with regard to jet breakup are that 
(1) pressure and density perturbations are a maximum at the wall and decrease to zero 
at the center, (2) tangential velocity perturbations are about 1. 5 times larger at the 
center than at the wall, and (3) radial velocity perturbations u r are zero at the wall but 
equal the u 0 perturbations at the center. The consequences of combining equal u 0 and 
u^ perturbations at the center is a rotating velocity vector of constant magnitude with no 
velocity perturbations. On approaching the wall, the steady component of this velocity 
vector decreases, and the perturbating component increases until it is entirely a pertur- 
bating component at the wall. This behavior is independent of the angular position. The 
effect of such acoustic properties on the jet dynamics can be shown by a parametric sur- 
vey of some of the variables. 

Effect of steady velocity difference and pressure amplitude at wall position. - At a 
wall position where velocity perturbations are in-phase and proportional to the pressure 
perturbations, the acoustic properties are similar to those of a traveling plane wave. 

The analytical and numerical results of the jet response at this position are shown in 
figure 10 for a pressure amplitude P' of 0. 1 and several values of the steady velocity 
difference ju^ _ u z | • The numerical results show the same periodic variation of the 
response factors with c of previously described for the analytical results. Quantita- 
tively, however, the results of the two analyses differ. 

The maximum peak value of (fig. 10(a)), obtained when |u^ - u z | is zero for 
both analysis, equals 1. 32 analytically and 1. 58 numerically. The separation between 
the analytical and numerical peaks with respect to c of also suggests a nonlinear velocity 
effect on jet breakup time. The numerical N T curves (fig. 10(b)) show an increase over 
the analytical results comparable to that for N^. For the constant P' conditions of 
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figure 10, the peak values decrease with an increase in | u ^ “ u z I and the difference be- 
tween the results of the two analyses decreases. 

The effect of both |u^ - u z | and P' on peak response is shown in figure 11. The 
peak values of N R near an cor of 3 jt are used for the comparison. When ju^ “ u z I is 
finite and constant, the response factor has an asymptotic minimum at low amplitudes. 

At high amplitudes, the analytical results approach a maximum asymptotically; however, 
the numerical result is more complex. As |u^ - u z | increases, pressure amplitude for 
the transition between these two asymptotic limits also increases. The minimum asymp- 
tote for the numerical results equals the analytical value of 1/y. The analytical value 
of 1 fy is the pressure sensitive response of the process when velocity oscillations are 
zero. The properties of the Nj response factor are similar to those N R characteristics. 

These results empahsize the importance of the in-phase velocity perturbation u R 
on the response of the jet atomization process. The magnitude of u R is proportional to 
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P\ At high amplitudes u R dominates |u^ - u z j and the response is significantly larger 
than for a dominating | | . In the case of the numerical results, emphasizing U R 
can almost double the response. A strong traveling wave or a low velocity difference 
environment, therefore, can cause a large dynamic response from the jet atomization 
process. 

The quantitative difference between the numerical and analytical results for a domi- 
nating in-phase velocity perturbation is caused by the nonlinear effects preserved in the 
numerical procedures. This nonlinear behavior is clearly evident in the cyclic variation 
in breakup time and atomization rate with time shown in figure 12. Two abrupt changes 



Figure 12. - Typical numerical solution for jet breakup time and 
atomization rate variations with time for pressure amplitude of 
0. 2. Steady velocity difference, 200 feet per second (61 m/sec); 
ratio of specific heats, 1.2; speed of sound, 5000 feet per sec- 
ond (1520 m/sec); maximum angular velocity, 452. 7 feet per 
second (138 m/sec). 


occur for each cycle of pressure oscillation. These are caused by the velocity pertur- 
bation. Both the positive and negative velocity peaks of the velocity perturbation enhance 
the breakup process. This dynamic behavior is characterized by and Np which 
express only the first harmonic content of the atomization rate perturbations. In fig- 
ure 12 the first harmonic content of the highly nonlinear atomization rate is obscure. In 
contrast, very small nonlinear effects were evident in the analytical atomization rate 
(fig. 7) where only the first harmonic content of the areodynamic forces on the jet were 
considered. Preserving the nonlinear content of these forces and eliminating other as- 
sumptions in the analytical procedures contribute to additional first-order content of the 
atomization rate perturbations. 

As previously shown in figure 9, velocity perturbations also affect the mean breakup 
time off. Figure 13 shows c of properties at the wall position in the traveling mode. 
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Figure 13. - Analytical normalized breakup time at wall position of traveling mode. 
Ratio of specific heats, 1.2 ; speed of sound, 5000 feet per second (1520 m/sec). 


As the pressure amplitude becomes large, substantial decreases in breakup time are 
predicted. The numerical solutions give comparable changes. Such changes in off 
cause significant changes in the response factors with amplitude because and Nj 
are periodically related to wr. 

This nonlinear effect of pressure amplitude on off introduces the possibility of non- 
linear instability limits for a rocket combustor. Referring to figure 8, a combustor with 
atomization limited combustion may be stable in an off region near 2rr radians, where 
Nj^ is negative. Velocity perturbation introduced by a bomb or other disturbance could 
reduce off to a region where Np is sufficiently high to overcome acoustic losses and 
cause instability. If the disturbance were sufficiently large to reduce off to near zero 
where N R is small, instability could occur during the decay of the disturbance. During 
decay (off increasing), the magnitude of increases. Therefore, if the peak re- 
sponse is larger than the acoustic losses, there should be some size of disturbance above 
which the combustor is always unstable. 

Effect of radial position. - The magnitude of the steady and perturbating content of 
the transverse velocity vector varies with radial position. The effect of such changes on 
breakup time cause the most significant changes in the response factors. Figure 14 
shows the change in breakup time with radial position. Breakup time (and jet length 
which is proportional to r) decrease on approaching the center. The fractional change 

decreases with an increase in lu 7 - u I . 

1 L z 1 
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Figure 14. - Analytical variation in normalized breakup 
time with radial position in traveling mode. Pressure 
amplitude, 0. 1 at radial position, 1.0; ratio of specific 
heats, 1.2; and speed of sound, 5000 feet per second 
(1520 'm/sec). 


Characteristic breakup 
time at wall, 



Radial position, r 

Figure 15. - Variation in in-phase response factor with 
radial position in traveling mode. Pressure amplitude, 
0. 1 at radial position, 1.0; and steady velocity differ- 
ence, 100 feet per second (30 m/sec). 
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It has been previously shown that the response factors and Nj vary periodi- 
cally with cot. Substantial changes in these response factors with radial position, there- 
fore, can be expected from changes in t with position. In fact, the response factors 
can exhibit periodic variations with radial position. Examples for several jet conditions 
are shown in figure 15. The value of N-^ can be positive at the wall and negative at the 
center, the converse, and many intermediate conditions. The values of Nj vary in a 
similar manner. Such changes in response with radial position imply that the manner 
in which jets are distributed within the cross-sectional area would significantly affect 
the average response. 

The dynamic behavior with radial position can deviate considerably from the ex- 
amples shown in figure 15. The variations of response with radial position, however, 
are contained within limiting values as indicated in figure 15. The limiting values at the 
wall are those previously discussed where low values of |u^ - u | increase the limiting 
value. At the center the limiting value is 1/y, the pressure sensitive value, for all 
conditions . 

The varied characteristics of the response factors with radial position preculdes 
any general characterization of the behavior. Only in the case of high axial velocity dif- 
ferences or low pressure amplitudes is the response relatively uniform throughout. 


Standing Transverse Acoustic Mode 


The first-order terms for the properties of the standing mode are given by 

P = p|l + AyJ 1 (o , )cos cut cos 

r jl + AJ-^c^cos cot cos 0] 


P g P gl 


Ji(a) 

u fl = Ac sin cut sin 9 

V a 


u r = Ac 


JJa) 

J o<“> - - 


a 


sin cut sin 9 


(38) 

(39) 

(40) 

(41) 


In the standing mode both u 0 and u r are out-of-phase with the pressure so that 

vl/2 


u i = ( u f + u ?) 
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(42) 



and the following identities are established with regard to equations (15) to (18) 


P' = AyJ 1 («)cos 9 




r 

u T = tan 9 { 2 + 
1 ya 


a 




H2 


Jj(«) 


- 2 



a = 1. 841 r 


( 43 ) 


In the standing mode, a pressure node exists along a diametrical path of 0 = n/2. Along 
the path of 9 = 0 the angular velocity is zero or at a nodal condition. As implied by 
these nodes, the standing mode differs from the traveling mode in that acoustic prop- 
erties vary with angular position as well as radial position. With regard to jet dynamics, 
however, complexity is somewhat reduced because velocity perturbations are restricted 
to out-of -phase perturbations that can only vary in amplitude between the limits 
±(c/«)(P ' Q=0 /y). 

The wall position at the angular velocity node (9 = 0) is unique with regard to jet dy- 
namics because velocity perturbations disappear and the density oscillations alone cause 
perturbations in the forces acting on the jet. This condition was studied in reference 14. 
In the analytical expressions for the response factors (eqs. (29) and (30)), the parameter 
V' is equal to one and is independent of the steady axial velocity difference and the pres- 
sure amplitude. The maximum values of Np and Nj, therefore, are equal to about 
1/y. 

Numerical and analytical solutions for Np are compared in figure 16. Nonlinear 
effects preserved in the numerical analysis appear to have a small effect on response in 
that agreement with the analytical solution is good even at high amplitudes. Similar 
agreement is obtained for Nj. Jet response to density oscillations alone, therefore, 
give a maximum value of about l/y for the response factors. The result would be iden- 
tical for a velocity node position in a longitudinal or any other mode. 

At any other position within the standing mode, the maximum values of N^ and Nj 
are affected by velocity perturbations restricted to out-of-phase properties. The effect 
of an out-of -phase velocity is shown in figure 17. Analytical and numerical solutions 
for Np are compared for |u^ - u z | = 0 at a quarter wave position (9 = tt/4) at the wall. 
At cot = 3jt radians, the maximum analytical response is l/2y (0.417) or 1/2 the pres- 
sure sensitive value at these dominating out-of-phase velocity conditions. The numeri- 
cal solution shows a maximum response of only 0. 338 at near 3 j r radians. This shows 
that an out-of-phase velocity depresses the response, and it was previously shown that 
an in-phase velocity enhanced the response. It can be concluded that a Uj velocity has 
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Figure 16. - Response at velocity node in standing mode. 


Type of solution 



Figure 17. - Response at angular position of 7r/4 radians and radial 
position of 1.0 in standing mode. Steady velxity difference, 0; 
pressure amplitude, 0. 1; speed of sound, 5000 feet per second 
(1520 m/sec); ratio of specific heats, 1.2. 



t 


a stabilizing effect and a velocity a destabilizing effect on a combustion system sen- 
sitive to jet atomization. 

The depressive effect of an out-of-phase velocity on response is more clearly shown 

in figure 18. Analytical values of (Nr, ) are shown as a function of P* and lu. - u I 

V R 4nax 1 1 z 1 

for the quarter wave wall position. The characteristics shown in figure 18 are an inver- 
sion of a similar variation for the traveling wave shown in figure 11. For a given value 
of | u £ - u z | in this standing mode has an asymptotic maximum at low pressure am- 
plitudes and an asymptotic minimum at high amplitudes. The suppressive effect of an 
out-of-phase velocity is most severe at high pressure amplitudes and low values of 
|u 7 - u | . Numerical solutions gave similar trends but the suppressive effect is more 
severe as shown in figure 18. 

Similar properties of (n r ) are exhibited at all other positions within the standing 

mode except for the constant value of (n r j of 1 / y occurring at the velocity node. 

'max 

The pressure node is no exception even though pressure oscillation disappears. In a sys- 
tem analysis for stability, however, the response at a pressure node usually becomes 
insignificant because system integrals essentially weight the response by the pressure 
amplitude. 

Generalizing the response in the standing mode is difficult. Response varies with 
both angular and radial position and causes complex patterns of response factors distri- 
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Figure 18. - Maximum in-phase response factor at angular position of tt/ 4 radians and radial position of 1.0 in 
standing mode. Speed of sound, 5000 feet per second U520 m/sec); ratio of specific heats, 1.2. Maximum 
values near u/t = 3tt used for comparison. 
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buted across the cross section. The maximum average value of either or Nj in the 
standing mode, however, will always be less than for the traveling mode. 


Radial Acoustic Mode 

The first order terms for the properties of the radial mode are given by 


P = p[l + AyJ Q (/3)cos cot] 

(44) 

P g = P g [i + AJ Q (/3)cos o>tJ 

(45) 

u r = -Ac/3jj(/3)sin cot 

(46) 


where 


/3 = 3.83 r 

Radial velocities are out-of-phase with pressure and there are no angular velocities. 
The identities become 


11 

c 

cP' 

p 1 

y 



-'i 


> 


(47) 


P’ = AyJ o (0) 


Maximum values of and Nj for the radial mode are contained within the same 

limits as for the standing mode. At the wall and center, where velocity perturbations 

vanish, (n^) is l/y. At other positions (Nr) can *> e reduced to about ^ y an- 
rricix max 

alytically and 0. 33 numerically because of the suppressive effect of an out-of-phase ve- 
locity component. 

As in the case of the other modes, velocity perturbations can cause breakup time to 
vary with radial position and substantially effect the values of and Nj. Figure 19 
compares changes in t with radial position for the different modes. The comparison is 
for | u £ - u z | equal to 100 feet per second (30. 5 m/sec) and a maximum pressure ampli- 
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Figure 19. - Comparison of normalized breakup time vari- 
ations for various acoustic modes. Maximum pressure 
amplitude, 0. 1 at radial position of 1.0 for all modes. 
Steady velocity difference, 100 feet per second (30 m/sec). 


tude at the wall of 0. 1 for all modes. For the radial mode off attains a minimum value 
at the pressure node. The traveling mode exhibits the smallest total change in c or when 
considering both radial and angular position. Considering this distinction in addition to 
the lower peak values of N-^ in the standing and radial modes further implicates the 
traveling mode as potentially most unstable. Large changes in off and the accompanying 
changes in can severely degrade the potential of the standing and radial modes for 
high response. This does not imply that space averaged response factors are always the 
largest in traveling modes. 


EFFECT OF DISTRIBUTED PROPERTIES 

The preceding analyses have shown the dynamic properties of a discrete breakup 
process in which all elements of a jet breakup when distortion attains a constant critical 
value. Breakup of an actual jet, however, rather than being discrete is distributed in 
both time and space as shown in figure 6. Such distributed properties will affect the dy- 
namic response. For example, figure 6 shows the time interval during which breakup 
occurs to be proportional to the mean breakup time. If the mean breakup time is large, 
the large time interval of breakup and mass release could encompass several cycles of 
acoustic oscillation. Response for such a condition would be small because the time 
correlation between mass release rate and pressure oscillations (a basic property of re- 
sponse factors) is small. Similarly, increasing the frequency would eventually decrease 
the correlation and reduce the response. This implies that the response factors Nj^ and 
Nj should decrease with increasing values of «r. The response of discrete breakup 
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(eq. (29) and (30) and fig. 7) does not exhibit such properties. Rather, the response 
factors beyond an wr of about 4 jt radians become simple periodic functions of con- 
stant amplitude. Dynamic properties of a discrete breakup process, therefore, do not 
characterize those of the actual process when cot is large. 

In the following analysis the modifying effects of a distributed process on the dy- 
namic properties of a discrete process are examined in an attempt to more realistically 
model the actual characteristics of jet breakup. 

For the purpose of analysis, a multiplicity of identical jets are assumed in a com- 
mon environment. The dynamic properties of each jet are those of a discrete breakup 
process. The critical value of distortion e identified with breakup, however, is a dis- 
tinguishing feature of each jet. The population density of jets as a function of the breakup 
time r is chosen to be the mass release distribution with r shown in figure 6. 

The mass release rate in-phase with the pressure oscillation for any jet is 
given by equation (29). Since pressure disturbances are identical for each jet in the 
common environment, the average value of the mass release rate in-phase with pressure 
for the multiplicity of jets is the integral average of the individual values of Nj^. Ap- 
proximating the population density of jets by a uniform distribution with outer boundaries 
of ±cr, the equation for the average response is given by 


jl n r j(ua) k 

V’ R 2 acor XrU-crt ' 


sin c or 


- COS COT )d(jOT 


o>T(l-a) 


COT 


(48) 


The solution to equation (48), where the integration of the first term is given by a 
series expression, is 


2— N = 
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When cr is small, the solution can be approximated by 


-r— Nr, = 


sin ocot /sin cor 
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CTCOT \ COT 


- COS COT 


Similarly, the average value of Nj is approximated by 


(50) 
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( 51 ) 


in t = si g^ T A - eos cor _ sin 

V T 1 (JCjOT \ COT , 

The term sin acor/crcor is a converging period function which modifies the response 
expressions for the discrete breakup process. A value of cr of 0. 15 encompasses 
50 percent of the mass release for the distribution curve of figure 5. Using a cr of 0. 15 
to characterize the distribution, the modified response function given by equations (50) 
and (51) are shown in figure 20. With the modified response functions, the peak values 



Figure 20. - Response for distributed jet breakup process. Deviation parameter, u)o, 0. 15. 


decrease with an increase in c of, the type of behavior expected from an actual jet atom- 
ization process. 

Accepting the modified response functions as more characteristic of the actual proc- 
ess does not particularly alter the previous discussion of the discrete process in specific 
acoustic modes. The initial maximum value of these periodic response functions is not 
significantly affected in the modified functions and is comparable to the value at oor = 3tt 
used previously in figures 11 and 18 as an index of the potential response. 

The properties of a distributed process emphasize the difficulties to be expected in 
experimental observations of jet dynamics. The breakup is obscure because it is distri- 
buted in both length and time. Distributed breakup also reduces the magnitude of any 
change in atomization rate. In addition, the dynamic properties of primary importance 
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are related to the first harmonic content of these distributed properties. Indeed, the 
prospect for reliable dynamic data from direct measurements of jet length properties 
seem remote. 


EFFECT OF VORTEX FLOW 

Vortex flow was a primary variable in the two-dimensional combustor studies of 
liquid jet behavior during instability shown in figures 2 to 4. In these studies, the degree 
of instability was related to the degree of vortex flow. This implies that dynamic re- 
sponse of the combustion process is affected by a steady angular flow. Analytical 
studies (refs. 15 to 17) have shown that a steady velocity can enhance the dynamic re- 
sponse of a velocity sensitive process. An analysis of the effect of a steady angular ve- 
locity on the jet breakup process will show that a significant increase in dynamic re- 
sponse can also be expected. 

Limiting the analysis to a wall position in a traveling transverse-acoustic mode 
(simulating the two-dimensional combustor studies), the forcing function f(t) for equa- 
tion (13) takes the form 


f (t ) = Pg(|u z - u z | 2 + u 0 2 ) 


(52) 


where 


P(t) = P(1 + P’ cos cat) 


P g (t) = P g jl + - P’ cos cat 


u„ (t) = u„ + — — cos cat 

0 0 Ol y 


Using these properties and extracting the first harmonic content of f(t) gives 


«« = ^rms)j( 1 + ^^co S cct) 


(53) 
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where 



( 54 ) 


(55) 


The two response factors are of the same functional form as in equations (29) and (30) 
Including the modifying effect of a distributed process, they are 

N _ s in go >t /sin lot _ 

" y OLOt \ lot 

XT v 0 sin <7 lot (l - cos lot 

Nj = — = 

y cor \ Poor 

The mean breakup time wr is given by 




COT ~ 



(58) 


The converging periodic nature of both and Nj with cot is identical to that 
previously discussed and shown in figure 20. The steady angular velocity u^, however, 
affects the modifying factor Vg/y and thereby significantly affects the magnitude of the 
response. Mean breakup time is also sensitive to the angular velocity. These affects 
are best illustrated by a specific example. 

The conditions citied in reference 8 for the oxygen jet breakup process shown in fig- 
ure 3 are 
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Jet diameter, D, in. ; m 

Average angular velocity, u fl , ft/sec; m/sec 

Liquid density, p 7 , lbm/ft ; kg/m 

Gas density, p , lbm/ft ; kg/m 

Liquid velocity, ju.^, in. /sec, m/sec . . . . 
Average radial velocity, in. /sec; m/sec . . 

Frequency, o>, rad/sec 

Speed of sound, c, ft/sec; m/sec 

Pressure ratio, AP / P 

Ir r 


0.028; 0.00071 
. . . 235; 72 
. . 75; 1200 

. . 0 . 1 ; 1 . 6 
. . 360; 9. 2 

. . 1140; 29 

.... 27 300 
. . 4720; 144 
.... 0. 667 


For u h u in equations (54) and (55), a characteristic time cot of 2860 radians and a 
r z 

jet length of 0. 72 inch (0. 018 m) are predicted for the case of no angular velocity, and 


e = 3. 5. Figure 21 shows the effect of a steady tangential velocity on breakup time and 


jet length for several pressure amplitudes. As the tangential velocity increases to 
300 feet per second (72 m/sec), the breakup time is reduced to the period of about 1^ 
cycles of acoustic oscillation where a significant response from the breakup process is 
anticipated. The peak magnitude of the response factors is approximated by Vj/y, 



0 50 100 150 _200 250 300 

Mean tangential velxity, Ug, ft/sec 



0 20 40 60 80 100 

Mean tangential velocity, Ug, ml sec 

Figure 21. - Calculated effect of tangential velxity and pressure am- 
plitude on jet length and mean breakup time for test conditions of 
reference 8. 
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0 20 40 60 80 

Mean tangential velocity, uq, m/sec 


Figure 22. - Calculated effect of tangential velocity and pressure am- 
plitude on response magnitude for test conditions of reference 8. 


which is also dependent on u^-. This is shown in figure 22. The amplifying effect of 
on response is substantial. The response factors can be increased by more than an 
order of magnitude over that for no steady tangential velocity when the response factor 
is near unity. 

The specific pressure amplitude for the jet conditions of this example are not di- 
rectly applicable to the dynamic analysis because the pressure waveform is strongly non- 
linear. A sinusoidal wave amplitude of P' = 0. 3, which is the approximate first harmonic 
content of the actual wave, may be representative. For P* <= 0. 3 and u„ = 235 feet per 

^ U 

second (72 m/sec), an cur/2 n of 1. 65 and V^/y of 10. 8 are given by figures 21 and 22. 
These values prescribe response factors of 4. 0 for N R and 6. 5 for Nj which would 
usually imply a highly unstable condition. 

The effects of u 0 on N R are more clearly shown in figure 23. For a constant 
amplitude of 0. 3, an increase in u^ is shown to cause periodic variations in N R of in- 
creasing amplitude. The major region of positive response occurs for tangential veloc- 
ities above 220 feet per second (67 m/sec) and encompasses the reference 8 test condi- 
tion. 

The effect of pressure amplitude on response for a tangential velocity of 235 feet 
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Mean tangential velocity, u Q , ft/sec 



0 20 40 .60 80 


Mean tangential velocity, tig, m/sec 

Figure 23. - Calculated effect of tangential velocity on in-phase re 
sponse factor for pressure amplitude of 0.3 and other jet condi- 
tions of reference 8. 



Pressure amplitude, P' 

Figure 24. - Calculated effect of pressure amplitude on in-phase re 
sponse factor for mean tangential velocity of 235 feet per second 
(71 m/sec) and other jet conditions of reference 8. 



per second (72 m/sec) is shown in figure 24. Positive response and potential instability 
is indicated for the reference 8 pressure amplitude of 0. 3. Figures 23 and 24, imply 
unstable combustion in regions other than the test condition of reference 8. Experimen- 
tally (refs. 7 and 14), pressure amplitude increased uniformly with u^. A stability 
analysis of the complete system is required, however, to establish the relation between 
equilibrium amplitude and combustion response. 


SUMMARY OF RESULTS 

The dynamic response of a liquid jet atomization process within traveling, standing, 
and radial transverse acoustic modes was derived by both theoretical and numerical anal- 
yses. The response was characterized by two response factors: the in-phase N R and 
out-of-phase Nj components of atomization rate perturbations with respect to pressure 
perturbations. A characteristic breakup time a >t (the product of jet breakup time and 
frequency) was an important correlating parameter for the response factors. Results of 
the analysis are as follows: 

1. For all cases, the response factors N R and Nj were periodic functions of the 
characteristic breakup time cot. The amplitude of the periodicity converges with an in- 
crease in cot when a breakup process distributed along the jet length is considered. 

The maximum values of N R and Nj occur at values of cot of about n and 3rr/2 ra- 
dians, respectively, with secondary peaks occurring at intervals of about 2ir radians. 
Convergence of the periodic function implies that large increases in jet diameter or 
oscillatory frequency would stabilize a combustion system. 

2. The magnitude of the peak values of the response factors varied with acoustic 
conditions. Values of N R as large as 1. 25 for the traveling mode and 0. 833 for the 
standing and radial mode were obtained by theoretical analysis. The numerical analysis, 
which preserved all nonlinear effects, gave a larger value of 1. 51 for the traveling mode. 
Comparable values were obtained for Nj. 

3. Variations in the peak values of the response factors are caused by the velocity 
sensitivity of the atomization process. Pressure sensitivity alone gave a peak value of 
0. 833. A velocity perturbation in-phase with the pressure elevated the peak value and is 
destabilizing to a combustion system. An out -of -phase velocity perturbation reduced the 
peak value and is stabilizing. A large steady relative velocity difference between the 
liquid and gas suppresses these velocity effects and causes the response to approach the 
pressure sensitive value. 

4. Spacial variations of acoustic properties within a given mode can cause large 
variations in local dynamic properties. The most severe effects on response factors are 
caused by the effect of acoustic properties on mean jet breakup time, which varies in- 
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versely with the magnitude of the local velocity perturbations. Response factors are 
periodic functions of breakup time and thus vary spatially throughout the mode. Such 
varied response precludes any simple characterization of the average response in a given 
mode. Potentially, however, the maximum average response for the traveling mode is 
larger than for either the standing or radial modes. Average response also depends on 
the uniformity of the injection pattern and may account for the effect of nonuniform in- 
jection on stability observed experimentally. 

5. The amplifying effect of a steady vortex velocity on the response of the jet atom- 
ization process is substantial. Peak response factors can be increased by an order of 
magnitude. The result may explain the instability encountered in a circular research 
combustor with vortex flow. Quantitative evaluations for one test condition showed jet 
response characteristics consistent with the experimental observations. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 9, 1969, 

128-31-51-02-22. 


40 



APPENDIX - SYMBOLS 


A amplitude coefficients, eqs. (37), 

(43), and (47), dimensionless 

a acceleration, ft/sec 2 ; m/sec 2 

c speed of sound, ft/sec; m/sec 

D jet diameter, ft; m 

J Bessel function of the first kind 

n 

of order n 

N response factor, component of 

w'/P’, dimensionless 
2 

P pressure, psi; N/m 

r fraction of combustor radius or ra- 
dial position, dimensionless 

t time, sec 

u velocity, ft/sec; m/sec 

V magnitude of relative velocity 

vector, ft/sec; m/sec 

w mass flow rate, lbm/sec; kg/sec 

a radius parameter, 1. 841 r 

P radius parameter, 3.83 r 

y ratio of specific heats, dimen- 

sionless 

6 jet distortion, ft; m 

e relative distortion, 6/D 

9 angular position, radians 


p density, lb/ft 2 ; kg/m 

a deviation from mean, dimen 

sionless 

t jet breakup time, sec 

o> frequency, rad/sec 

Subscripts: 
g gas 

I imaginary or out-of -phase 

component 

l liquid 

p-p peak- to -peak 

R real or in-phase component 

r radial component 

rms root mean square 

ss steady-state 

z axial component 

9 angular component 

0 initial 

Superscripts 

~ mean or average value 

maximum value 

* perturbation about mean 
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